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We describe the structure of the extended Clifford Group (defined to be the 
group consisting of all operators, unitary and anti-unitary, which normalize the 
generalized Pauh group (or Weyl-Heisenberg group as it is often called)). We 
also obtain a number of results concerning the structure of the Clifford Group 
proper (i.e. the group consisting just of the unitary operators which normalize 
the generalized Pauli group). We then investigate the action of the extended 
Clifford group operators on symmetric informationally complete POVMs (or 
SIC-POVMs) covariant relative to the action of the generalized Pauli group. 
We show that each of the fiducial vectors which has been constructed so far 
(including all the vectors constructed numerically by Renes et al) is an eigen- 
vector of one of a special class of order 3 Clifford unitaries. This suggests a 
strengthening of a conjecture of Zauner's. We give a complete characterization 
of the orbits and stability groups in dimensions 2-7. Finally, we show that 
the problem of constructing fiducial vectors may be expected to simplify in the 
infinite sequence of dimensions 7, 13, 19, 21, 31, ... . We illustrate this point by 
constructing exact expressions for fiducial vectors in dimensions 7 and 19. 
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1. Introduction 

The statistics of an arbitrary quantum measurement are described by a positive 
operator valued measure, or POVM (Davies 1 , Busch et al 2j, Peres [31, Nielsen 
and Chuang 4 and references cited therein) . Suppose the measurement has only a 
finite number of distinct outcomes. Then the corresponding POVM assigns to each 
outcome i the positive operator Ei with the property that Tv{Eip) is the probability 
of obtaining outcome i (where p is the density operator). Since Ti{Eip) — 1 for 
all p we must have J^iJ^i — 1- 

A POVM is said to be informationally complete if the probabilities Tv{Eip) 
uniquely determine the density operator p. The concept of informational complete- 
ness is originally due to Prugovecki |S] (also see Busch 0, Busch et al |2J, d'Ariano 
et al 0, Flammia et al [HI, Finkelstein [H], and references cited therein). It has an 
obvious relevance to the problem of quantum state determination. It also plays an 
important role in Caves et o/'s ^HEIEIEI Bayesian approach to the interpre- 
tation of quantum mechanics, and in Hardy's |14l I15| proposed axiomatization. 

Suppose the Hilbert space has finite dimension d. Then it is easily seen that an 
informationally complete POVM must contain at least (P distinct operators Ei. An 
informationally complete POVM is said to be symmetric informationally complete 
(or SIC) if it contains exactly this minimal number of distinct operators and if, in 
addition, 

(1) \Ei is a one dimensional projector for all i and some fixed constant A. 

(2) The overlap Tr{EiEj) is the same for every pair of distinct labels i,j. 

It is straightforward to show that this is equivalent to the requirement that, for 
each i, 

E^ = \ m (V'd (1) 

where the (P vectors \'4>i) satisfy 

um\ = \^^ '.y. (2) 

SIC-POVMs were introduced in a dissertation by Zauner ^^Ij and in Renes et 
al ^Tj. Wootters ^1], Bengtsson and Ericsson ^l|5n| and Grassl |2] have made 
further contributions. There appear to be some intimate connections with the 
theory of mutually unbiased bases ^11221 ESI j finite affine planes ^JEIED], and 
polytopes [mill. 

If SIC-POVMs existed in every finite dimension (or, failing that, in a sufficiently 
large set of finite dimensions) they would constitute a naturally distinguished class 
of POVMs which might be expected to have many interesting applications to quan- 
tum tomography, cryptography and information theory generally. They would also 
be obvious candidates for the "fiducial" or "standard" POVMs featuring in the 
work of Fuchs [TS| and Hardy [11 US]. 

The question consequently arises: is it in fact true that SIC-POVMs exist in every 
finite dimension? The answer to this question is still unknown. Analytic solutions 
to Eqs. l|21l have been constructed in dimensions 2, 3, 4, 5, 6 and 8. Moreover Renes 
et al jl7j have constructed numerical solutions in dimensions 5 to 45 (the actual 
vectors can be downloaded from their website ^j). So one may plausibly speculate 
that SIC-POVMs exist in every finite dimension. But it has not been proved. 

The SIC-POVMs which have so far been explicitly""^ described in the literature are 
all covariant under the action of the generalized Pauli group (or Weyl-Heisenberg 



"'^Renes et al 1171 mention that they have constructed numerical solutions which are covariant 
under the action of other groups, but they do not give any details. 
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group, as it is often called). It is therefore natural to investigate their behaviour 
under the action of the extended Clifford group. The Clifford group proper is 
defined to be the normalizer of the generalized Pauli group, considered as a subgroup 
of U(d) (the group consisting of all unitary operators in dimension d). It is relevant 
to a number of areas of quantum information theory, and it has been extensively 
discussed in the literature [23 EH EZl EH EH EOl EB ■ Its relevance to the SIC- 
POVM problem has been stressed by Grassl As Grassl notes, it is related to 
the Jacobi group ^22,, which has attracted some notice in the pure mathematical 
literature. We define the extended Clifford group to be the group which results 
when the Clifford group is enlarged, so as to include all anti-unitary operators 
which normalize the generalized Pauli group. As we will sec, this enlargement is 
essential if one wants to achieve a full understanding of the SIC-POVM problem. 

In Sections Elfil we give a self-contained account of the structure of the extended 
Clifford group. In the course of this discussion we obtain a number of results 
concerning the structure of the Clifford group proper which, to the best of our 
knowledge, have not previously appeared in the literature and which may be of 
some independent interest. 

In Sectional we define and establish some of the properties of a function we call 
the Clifford trace. We also identify a distinguished class of order 3 Clifford unitaries 
for which the Clifford trace = — 1. We refer to these as canonical order 3 unitaries. 

In Section El we analyze the vectors constructed numerically by Renes et al ^7] 
(RBSC in the sequel) in dimension 5-45. We show that each of them is an eigen- 
vector of a canonical order 3 Clifford unitary. This suggests the conjecture, that 
every CP fiducial vector is an eigenvector of a canonical order 3 unitary. We also 
show that, with one exception, the stability group of each RBSC vector is order 3 
(the exception being dimension 7, where the stability group is order 6). 

In Section [7| we show that RBSC's results also support a strengthened version 
of a conjecture of Zauner's TB' (also see Grassl j21j). 

In SectionlHlwe use RBSC's numerical data, regarding the total number of fiducial 
vectors in dimensions 2-7, to give a complete characterization of the orbits and 
stability groups in dimensions 2-7. Our results show that in each of these dimensions 
every fiducial vector covariant under the action of the generalized Pauli group is 
an eigenvector of a canonical order 3 Clifford unitary. We also identify the total 
number of distinct orbits. It was already known ^| |2] that there are infinitely 
many orbits in dimension 3, and one orbit in dimensions 2 and 6. We show that 
there is, likewise, only one orbit in dimensions 4 and 5, but two distinct orbits 
in dimension 7. We also construct exact expressions for two fiducial vectors in 
dimension 7 (one on each of the two distinct orbits) . 

RBSC's numerical data may suggest that, after dimension 7, the stability group 
of every fiducial vector has order 3. In SectionlHlwe show that there is at least one 
exception to that putative rule by constructing an exact expression for a fiducial 
vector in dimension 19 for which the stability group has order > 18. 

Our construction of exact solutions in dimensions 7 and 19 was facilitated by 
the fact that in these dimensions there exist canonical order 3 unitaries having a 
particularly simple form. In Section [Till we show that a similar simplification occurs 
in every dimension d for which (a) d has at least one prime factor — 1 (mod 3), (b) 
d has no prime factors = 2 (mod 3) and (c) d is not divisible by 9. In other words, 
it happens when d ^ 7, 13, 19, 21, 31, .... 

2. Fiducial Vectors for the Generalized Pauli Group 

The SIC-POVMs which have been constructed to date all have a certain group 
covariance property. Let G be a finite group having d"^ elements, and suppose we 
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have an injective map g ^ Ug which associates to each g G G a unitary operator 
Ug acting on d-dimensional Hilbert space. Suppose that for all g, g' 

UgUg,=e^i^^'Ugg, (3) 

where e'^ss' is a phase (so the map defines a group homomorphism of G into the 
quotient group U((i)/ Uc(d), where Uc(fi) is the centre of U(d)). Finally (and this, 
of course, is the difficult part) suppose we can find a vector {ip) S such that 
{tp\tp) = 1 and 

|(V'|J7,|V)| = ^=^ (4) 
for all 5 7^ e (e being the identity of G). Then the assignment 

Eg^^Ug\^P){ml 

defines a SIC-POVM on C^. The vector is said to be a fiducial vector. 

To date attention has been largely focussed on the case G = (Z^)^, where 
is the set of integers 0, 1, . . . , d — 1 under addition modulo d (although there is 
numerical evidence that fiducial vectors exist for other choices of group ^Tj). That 
is also the case on which we will focus here. 

To construct a suitable map (Z^)^ U((i), let |eo), |ei), . . . |ed_i) be an or- 
thonormal basis for C^, and let T be the operator defined by 

f|eO =a;ne.) (6) 
where co — e^'^^l'^ . Let S be the shift operator 

A, \ /|er+i) r = 0,l,...,d-2 

Then define, for each pair of integers p = {v\iVi) G Z^, 

bp = rP^P^SP'fP'' (8) 

where t — —e^'-l'^ (the minus sign means that r'' = 1 for all d, thereby simplifying 
some of the formulae needed in the sequel). We have, for all p, q G 1? , 

b; = b^p (9) 

bpb^ = r<P-i>^p+q (10) 

and 

n _ J^P if d is odd 

^P+rfq- |(_^)(p,q)^^ if d is even ^ ' 

where (p, q) is the symplectic form 

(P,q) =P29i -Pig2 (12) 

Consequently the map p E {'^df' bp e U(d) has all the required properties. 
The operators Dp are sometimes called generalized Pauli matrices. So we will say 
that a vector e C'' is a generalized Pauh fiducial vector, or GP fiducial vector 
for short, if it is a fiducial vector relative to the action of these operators: i.e. if 

(■(/'IV-') = 1 and 

Mbpm = ^^ (13) 



Vd + l 



for every p e 7^ (mod d) . 
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The set of operators Dp is not a group. However, it becomes a group if we 
allow each Dp to be multiplied by an arbitrary phase. We will refer to the group 
GP{d) = {e^^Dp : ^ e M, p G Z^} so obtained as the generalized Pauli group^. 

We now want to investigate the normalizer of GP(d): i.e. the group G{d) con- 
sisting of all unitary operators U € U(d) with the property 

U GP {d)U^ ^GP{d) (14) 

The significance of this group for us is that it generates automorphisms of GP(d) 
according to the prescription 

P UPij^ (15) 

Consequently, if j-i/') is a GP fiducial vector, then so is tJ\'ip) for every U S G{d). 

The group G{d) is known as the Clifford group, and has been extensively dis- 
cussed in the literature [23 EEl EH EHl ES EDI EH ■ Its relevance to the SIC-POVM 
problem has been stressed by Grassl However, none of these accounts derive 
all the results needed for our analysis of the RBSC vectors. In the interests of 
readability we give a unified treatment in the next section. 

3. The Clifford Group: Structure, and Calculation of the Unitaries 

We begin with some definitions. Let 

— Id if d is odd 

I 2a II ct IS even 

Let SL(2,Z^) be the group consisting of all 2 x 2 matrices 

such that a, /3, 7, 5 €E and a5 — (3^ — 1 (mod d). Note that inverses exist in this 
group because the condition a6 — (3^ — 1 (mod d) implies 

.7 S) [-J a J - [O 1, 

in arithmetic modulo d. 
We then have 



(18) 



Lemma 1. For each unitary operator U £ G{d) there exists a matrix F G SL(2. Z^) 
and a vector x S (Z^)^ such that 

UDpU^ ^J^^^P^Dpp (19) 
for all p G Z^ (where uj — — e^^'/'^, as before). 

Proof. If [/ G G{d) it is immediate that there exist functions / and g such that 

UDpU^ = e''(P)^/(p) (20) 
for all p G Z2. It follows from Eq. ^ that 

for all p, q G Z^. Consequently 



^Also known as the Weyl-Heisenberg group. Our definition is, perhaps, slightly unconventional. 
It would be more usual to define GP((i) = {r"Dp: n G Z,p £ 7p} — i.e. the subgroup generated 
by the operators Dp. 
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which imphes /(p + q) = /(p) + /(q) (mod d). We may therefore write 

/(p) = F'p + dh{p) (23) 

for some matrix F' and function h. Inserting this expression in Eq. (|2()(l gives, in 
view of Eq. I^IJ , 

UD,U^ . e-(P)i)^,p+,,(p) = 1'"!'!''^;^, , ^ "'^^ (24) 

With the appropriate definition of g' this means 

UDpU'' = e'a'{v)f)p,^ (25) 
for all p. Repeating the argument which led to Eq. H22(l we find 

gig'(p+q)-3'(p)-g'(q)^(p,q>-(F'p,_F'q) _ ^ (26) 

Interchanging p and q gives 

gig'(p+q)-g'(p)-g'(q)^-(p,q> + (F'p,F'q> _ ^ (27) 

We consequently require 

^(p,q>-(F'p,F'q> ^ .^2((p,q)-(F'p,F'q)) ^ ^ (33) 

for all p, q. It is readily verified that {F'p,F'q) ~ (Det i^')(p, q). We must there- 
fore have 

Det F' = 1 (mod d) (29) 
If d is odd, or if d is even and Deti^' = 1 (mod d), we can find a matrix F € 
SL(2,Z^) such that F ^ F' (mod d). It then follows from Eq. ^ that Dpp = 
Dp'p for all p. 

Suppose, on the other hand, d is even and DetF' ^ 1 (mod d). Then Deti^' = 
d+1 (mod d). Write 



k1 ^ / 

We know aS — (3j — Dot F' is odd. So either a, 5 are both odd, or else f3,j are 
both odd. If a, (5 are both odd let 

A- / 

,0 



A-fi n) (31) 



while if /?, 7 are both odd let 



A=(2 J) (32) 



Then Det(i^' + dA) = 1 (mod d). We can therefore choose a matrix F e SL(2, Z^) 
such that F ^ F' + dA (mod d). Inserting this expression in Eq. H25(l we have, in 
view of Eq. l(TT|l . 

UDpU^ = e^^'(P)^(F_rfA)p = e'^'(P)(-l)^^P'''P^^Fp (33) 

We conclude that there is, in every case, a function g" and a matrix F £ SL(2, Z^) 
such that 

;7£>pC/^ e'3"(P^DFp (34) 

for all p. 

It remains to establish the form of the function g". We note, first of all, that it 
follows from Eqs. ^ and (^HJ that 

{DpY = Ddp = t'^'p^p^S'^p^T'^p^ = 1 (35) 



for all p (because S'^ ^J"^ = t'^'^ = !)• Consequently 

1 = U{bpyU^ = {ubpU^y = e''^a"(p) [bp^Y = e"^9"(P) (36) 

for all p. We must therefore have e*^ '-P-' = w^(p) for some function g taking values 
in Zrf. Repeating the argument which led to Eq. H26() we find 

^9(p+q)-s(p)-9(q)^(p,q)-(f'p.-F'q) ^ i (-37^ 

We have (p, q) — (i^p, Fq) = (1 — Deti^)(p,q) — (mod d). Consequently 
^(p,q>-(Fp,_Fq> ^ I (because = 1) and so 

ff(p + q) - .9(p) + g(q) (mod d) (38) 

for all p, q. This implies ^(p) = (x'l P) (mod d) for for all p some fixed 'x^ E (Zd)^- 
Setting X = Fx' i and using the fact that {F^^x, P) = {X: Fp) (mod d) we conclude 

UbpU^ ^uj'^^'^pWfp (39) 
for all p. □ 

We now want to prove the converse of Lemma ^ That is, we want to prove 
that, for each pair F G SL(2,Z^) and x G (^d)^ there is a corresponding operator 
U e C{d). We also want to derive an explicit expression for the operator U (this 
has, in effect, already been done by Hostens et al however, the formulae we 
derive are different, and better adapted to the questions addressed in this paper). 

We begin by focussing on a special class of matrices F. Let [ni, ri2, • ■ • , fir] denote 
the CCD (greatest common divisor) of the integers ni,n2, . . . ,nr. We define the 
class of prime matrices to be the set of all matrices 

G SL(2,Z^) such that (3 is non-zero and [/?, d] — 1 (so that (3 has a multiplicative 
inverse in Zj). We then have 



Lemma 2. Let 

7 5 



(41) 



he a prime matrix G SL(2, Z^). Let 



Vi 



F 



1 ^/3-(a.= -2..+5r^)|g^^^g^| (42) 



^/d 



r,s— 



(where l3 ^ G Zq is such that (3 ^(3 = 1 {mod d)). Then Vp is a unitary operator 
G C(d) such that 

VpbpVl = bp^ (43) 

for all p. 
Proof. Let 



and define 



S' = Df^^,^^ and T' = Df^p,s) (44) 
l/o) = 7^E(^T|eo) (45) 

^ r— 

It follows from Eq. that (T')"^ ^ 1. Consequently 

r'l/o) = l/o) (46) 
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It follows from Eq. (|10() that f'S' = uS'T'. So we can obtain a complete set of 
eigenvectors by laddering. Specifically, let 

\fr) = (^')''|/o) (47) 

for r = 1. Then 

f'\fr)=^^\fr) (48) 

for all r. Since 

(S") = 1 (as follows from Eq. jSHJ) we also have 

^'l/.> = l/.e.i) (49) 

for all r (where (Bd signifies addition modulo d). 

We next show that the vectors \fr) are orthonormal. It follows from Eqs. 
0, © and (Cnil that 

(f Oleo) = Dirp^rS)\eo) = rf"'^" Sf^'leo) (50) 

and consequently 

l/o> - (t^E^''^'^"'^) |eo) = (^Er^"'^'''^^'^^'^) |eo) (51) 

(where we have used the fact that t'* — 1). We need to be careful at this point, due 
to the fact that congruence modulo d need not imply congruence modulo d. Let qr 
be the quotient of /3r on division by d, and let tr be the remainder. So /3r = qrd + tr 
and 



r=0 



We have 



and 



gq^d+U ^ gt^ (53) 



(because r^"* = t''^ = !)■ Consequently 

The fact that [(3,d] = 1 implies that [P,d] = 1. It follows that, as r runs over 
the integers 0,1, . . . ,d ~ 1, so does tr (though not necessarily in the same order). 
Consequently 

l/") = ^E-'"'*'l^*) (56) 

It follows that 

ifrlfr) = (MiST'^iSrUo) - (/Ol/O) = 1 (57) 

The fact that (/r|/s) = when r 7^ s is an immediate consequence of the fact that 
\fr), \fs) are eigenvectors of T" corresponding to different eigenvalues. We conclude 
that 

{fr\fs)=5rs (58) 

as claimed. 
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We now want to calculate an explicit formula for \fr) when r > 0. It follows 
from previous results that 

\fr) = Drc.rMo) = — V r'^" ' (^)™ | (59) 

By an argument similar to the one leading to Eq. 1)56(1 we deduce 

= ^yr^"(^*'-2'-*+"'-')|e,) (61) 



(since — /37 = 1 (mod d)). Comparing with Eq. (|42fl we see that 

VF^Y.\fr){er\ (62) 
r=0 

which shows that Vp is unitary. Moreover, 

VfT1/;|/,) = VFTle.) - c^n/.) (63) 

for all r. Comparing with Eq. ^ we deduce VpTV^ = f'. Similarly VpSVl = S' . 
Hence 

VpDpVl = TPif^V>5'PirP=V'j; (64) 

^^(l-/37+a5)piP2^^p (66) 

- Dfp (67) 
for all p. □ 

To extend this result to the case of an arbitary matrix e SL(2,Z^) we need the 
following decomposition lemma, which states that every non-prime matrix can be 
written as the product of two prime matrices: 



F-iZ (68) 



Lemma 3. Let 

w- I 

be a non-prime matrix e SL(2, Z^) . Then there exists an integer x such that 5 + x[3 
is non-zero and [6 + x/3, d] ~ 1. Let x be any integer having that property, and let 

--(; <-) 

Then Fi, F2 are prime matrices G SL(2,Z^) such that 

F = F1F2 (71) 

Proof. Suppose, to begin with, that /3, 6 are both non-zero. Let k = [/3, 6]. We then 
have 

/3 = k(3o (72) 
d = kSo (73) 
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where [/?o,<5o] — 1- We also have [k,d] — 1 (because aS — = 1 (mod d)). The 
fact that (3q, Sq are relatively prime means we can use Dirichlet's theorem (see, for 
example, Nathanson ^| or Rose [Hi]) to deduce that the sequence 

Sq, (So + M, iSo + 2l3o),... (74) 

contains infinitely many primes. Consequently, there exists an integer x such that 
6o + x(3o and [^o + x/3o, d] = 1. The fact that k ^ and [k, d] = 1 then implies 
that S + x(3 and [S + xf3,d] = 1. The claim is now immediate. 

It remains to consider the case when /3, S are not both non-zero. If <5 = the 
fact that det F = 1 (mod d) would imply that /3 ^ and [(3, d\ = 1 — contrary to 
the assumption that the matrix F is non-prime. Suppose, on the other hand, that 
P = 0. Then the fact that detF = 1 (mod d) imphes that S ^0 and [S,d] = 1. So 
the claim is true for every choice of x. □ 

We can now deduce the following converse of Lemma ^ 

Lemma 4. Let {F,x) be any pair G SL(2, Z^) x (Z^)^. If F is a prime matrix 
define 

U = D^Vf (75) 

(where Vf is the operator defined by Eq. i4^ )- U ^ i^ non-prime choose two prime 
matrices Fi, F2 such that F = F1F2 (the existence of such matrices being guaranteed 
by Lemma\^ , and define 

U ^ D^Vf^Vf^ (76) 
(where Vp^ , Vp^ ^''^ operators defined by Eq. ^^). Then 

UDpU^ ^J^'^p'^Dfp (77) 

for all \i el? 

Proof. The claim is an immediate consequence of Eqs. ©, (|1U|I and Lemma |2 □ 

If [/, U' differ by a phase, so that U' = e*®[/, they have the same action on the 
generalized Pauli group: 

UbpU^ ^ U'bpU'^ (78) 
for all p. So the object of real interest is not the Clifford group itself, but the group 
C(d) / 1(d) which results when the phases are factored out. Here l{d) is the subgroup 
consisting of all operators of the form e'^/, where / is the identity operator and 
6* G R. The elements of C(c?)/I(d) are often called Clifford operations. 

Let SL(2, Z^) K (Zrf)^ be the semi-direct product of SL(2, Z^) and (Z^)^: i.e. the 
group which results when the set SL(2, Z^) x (Z^)^ is equipped with the composition 
rule 

{.Fi,Xi) ° (^2,X2) = {.F^F2,Xi+F^X2) (79) 
Then we have the following structure theorem, which states that G{d)/l{d) is nat- 
urally isomorphic to SL(2,Z^) k (Z^)^ when d is odd, and naturally isomorphic to 
a quotient group of SL(2,Z^) k (Z^)^ when d is even: 

Theorem 1. There exists a unique surjective homomorphism 

f : SL(2, Z^) K (Zrf)2 ^ C(d)/ 1(d) (80) 

with the property UIDpU^ = w<>'^-^P>i)pp for each U E f{F,x) and all p G Z^. 

// d is odd f is an isomorphism. If d is even the kernel of f is the subgroup 
Kf C SL(2, Z^) K (Zd)^ consisting of the 8 elements of the form 

'l-frd sd \ fsd/2Y 
td l + rd)\td/2)^ 

where r, s.t — or 1. 



(81) 
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Proof. An operator U G G{d) has the property 

UDpU^ = Dp (82) 

for all p if and only if it is a multiple of the identity. So it follows from results 
already proved that there is exactly one surjective map 

/ : SL(2, Z^) K (Zrf)2 ^ C(rf)/ 1(d) (83) 

such that UDpU'f = uj<^^^p>Dfp for each U € fiF,x) and aU p e Z^. The 
fact that / is actually a homomorphism is then an immediate consequence of the 
definitions. 

Let Kf be the kernel of /. Then (F, x) G Kf if and only if 

co<^^Fp>Dpp = Dp (84) 

for all p. For that to be true we must have F — 1 (mod d). If d is odd this 
implies Dpp = Dp for all p. Eq. H84|l then becomes ti;<^'P> = 1 for all p, implying 

Suppose, on the other hand, that d is even. The condition F = I (mod d) then 
implies that -F = 1 + dA, where A is a matrix of the form 



So the kernel is trivial, and / is an isomorphism as claimed. 



with ri,r2,s,t = or 1. Inserting this expression in Eq. (|84|l wc find, in view of 
Eqs. (IMni), that (F, x) e Kf if and only if 

1 = UJ<^'^P>DFpD^p = ^<X,P>^d<P,Ap> (gg) 

for all p. After re-arranging the condition becomes 

for all p. This is true if and only if ri = r2, Xi = sd/2 and X2 — td/2. □ 

We conclude with a result concerning the order of the group C{d) / l{d) which will 
be needed later on. Let ly^n, d) be the number of distinct ordered pairs {x, y) G (Z^)^ 
such that xy = n (mod d) . We then have 

Lemma 5. The order of the group C{d)/l{d) is 

|C(d)/I(d)| ^d^ [ ^zy(n,d)zy(n + l,d) J 

If d is a prime number this reduces to 

\C{d)/l{d)\ = d^{d^ - 1) 

Proof. We begin by showing that C{d)/l{d) and SL(2, Z^) k (Z^)^ have the same 
cardinality when considered as sets. This is true for all d, notwithstanding the fact 
that when d is even C{d)/l{d) and SL(2,Zd) k (Z^)^ are not naturally isomorphic 
as groups. 

The statement is immediate when d is odd. Suppose, on the other hand, that d 
is even. Let g : SL(2,Z2(i) — > SL(2,Zd) be the natural homomorphism defined by 

where [x]d denotes the residue class of x modulo d. It is easily seen that g is 
surjective. In fact, consider arbitrary 

F^(^ ^)gSL(2,Z,) (91) 
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Then aS — (3^ — 1 + nd for some integer n. If n is even then F G SL(2,Z^) and 
F = g[F). Suppose, on the other hand, that n is odd. Then either a or /? is odd. 
If a is odd F ~ g{F') where 



while if (3 is odd F ~ g{F") where 



(92) 



f)eSL(2,Z^) (93) 



^7 + d S 

Now let Kg be the kernel of 5. A matrix F E Kg if and only if 

\ td l + r2dj ^ ' 

where ri, r2, s, i = or 1 and (1 + rid){l + r2d) — std^ = 1 (mod 2c?). We have 

{I + rid){l + r2d) - std^ = I + {ri + r2)d (mod 2d) (95) 

(bearing in mind that d is even, so c?^ = (mod 2d)). Wc therefore require ri — r2. 
It follows that Kg consists of the 8 matrices of the form 



1 + rd sd 
td 1 + rd 



(96) 



where r, s,t = or 1. The fact that g is surjective and \Kg\ = 8 implies | SL(2, 7i2d)\ = 
8 1 SL(2, Zrf)|. In view of Theorem this means 

\C{d)/l{d)\ - i|SL(2,Z2rf)||(Zd)2| = |SL(2,Zd)||(Zd)2| = |SL(2,Zd) x {Zd)^\ 

(97) 

as claimed. 

We have shown that \C{d)/l{d)\ = |SL(2,Zd) k {Zd)^\ = d2|SL(2,Zd)| for all d, 
odd or even. It remains to calculate |SL(2, Z^) | . For each n G Z^ let M„ C SL(2, Z^) 
be the set of matrices 



a 
7 S 



(98) 



for which aS = n + 1 (mod d) and /37 = n (mod d). Clearly SL(2, Z^) U„=o 
and |M„| = i'(n, d)v{n + 1, d). It follows that 



|SL(2,Zd)| = ^j/(n,d)zy(n + l,d) (99) 



n=0 

Eq. H88(l is now immediate. 
If d is a prime number 

, „ f2d-l ifn==0(modd) , , 

v{n,d) = { / > (100) 

I a — 1 otherwise 

implying 

d-l 

^v{n,d)v{n + l,d) ^ d{d^ -I) (101) 

n=0 

Eq. H89|l is now immediate. □ 
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4. The Extended Clifford Group 

It can be seen from Eqs. (I^HHl and lfO|l that, if IV") = X)r=o V'rlcr) is a GP fiducial 
vector, then so is the vector — J2t=o ''^rl^r) obtained by complex conjugation. 
So to make the analysis complete we need to consider automorphisms of GP{d) 
which are generated by anti-unitary operators. 

An anti-linear operator is a map L: ^ with the property 

L {a\<j>) + (3\i^)) ^ a* L\(l)) + (3* L\^) (102) 

for all 10), G and all a,/3 G C. The adjoint is defined to be the unique 
anti-linear operator with the property 

(0|Lt|^) = {^\L\c^) (103) 

for all 10), lip) E C*. An operator U is said to be anti-unitary if it is anti-linear and 
C/tj7 = 1 (or, equivalently, UU^ = 1). 

We now define the extended Clifford Group to be the group EC((i) consisting of 
all unitary or anti-unitary operators U having the property 

U GP {d)U^ = GP{d) (104) 

Let us also define ESL(2,Z^) to be the group consisting of all 2 x 2 matrices 

such that a, /3,j,5 E l-^ aS — (3^ — ±1 (mod d). In the last section we showed 
that there is a natural homomorphism /: SL(2,Z^) k (Z^)^ G{d)/l{d). We 
are going to show that this extends to a natural homomorphism /e : ESL(2, Zq) k 
(Zd)^ ^EC{d)/l{d). 

Let J be the anti-linear operator which replaces components in the standard 
basis with their complex conjugates: 

i-=0 r=0 

Clearly ,V — J and ,P J = — 1. So J is an anti-unitary operator. Furthermore, 
it follows from Eqs. (jSHEJ that 

jbpP = (107) 



for all p, where 



J -{I 'a (108) 



So J e EC(d). Note that Det J = -1 (mod d), so J E ESL(2, Z^). 

Now let C*{d) be the set of anti-unitary operators E EC(c?) (so EC((i) is the 
disjoint union EC(d) = G{d) U C (d)). The mapping U t-^ JU defines a bijective 
correspondence between C*(d) and C{d). We can use this to prove the following 
extension of Theorem ^ 

Theorem 2. There is a unique surjective homomorphism 

/e : ESL(2, Z^) k {Zdf ^ EC(d)/ 1(d) (109) 
such that, for each {F, x) e ESL(2, Z^) k {Zdf and U E fE{F, x), 

ubpU^ ^Lo<^^^p>bFp (110) 

for all p. U is unitary if Dot F = 1 (modd) and anti-unitary if Det F = —1 (modd). 
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/e extends the homomorphism f defined in Theorem^ and has the same kernel. 
So /e is an isomorphism if d is odd, while if d is even its kernel is the subgroup Kf 
defined in Theorem^ 

Proof. Let tj be an arbitrary anti- unitary operator e C*(d). The fact that J, U are 
both anti-unitary means that JU is unitary. So JU S C(d). It then foUows from 
Theorem n that there exists € SL(2,Z^) x (Z^)^ such that 

{.JU)Dp{.JUy = u;<^''^'P>DF,p (111) 

for aU p. Define F — JF' and x — Jx' ■ In view of Eq. p()7|l . and the fact that 
= 1, we deduce 

UDpU^ = J{JU)Dp{JU)U^ = oj-<^'^^''P>Djp,.p = uj<^'^p>Dfp (112) 

for all p (where we have used the fact that < ^, J7 — < J^, Jr] > for all ^, ry). 
We have Det(F) = (Det J)(Det F') = -1, so {F,x) e ESL(2,Z^) k {Zdf. 

Reversing the argument we deduce the converse proposition: for each (F, x) G 
ESL(2,Z^) X (Zd)2, there exists U e EC(d) such that UDpW = uj<^'^p>Dfp for 
all p. The fact that an operator commutes with Dp for all p if and only if it is a 
multiple of the identity means that U is unique up to a phase. 

This establishes the existence and uniqueness of the homomorphism /e- The 
proof of the remaining statements is straightforward, and is left to the reader. □ 

Finally, we have the following result which, together with Lemma |5| enables us 
to calculate the order of FiC{d)/ l{d): 

Lemma 6. 

|EC(d)/I(d)| = 2|C(d)/I(d)| (113) 

for all d. 
Proof. The map 

Ul{d) .JUl{d) (114) 

defines a bijective correspondence between C*{d)/l{d) and G{d)/l(d). So the set 
C {d)/l{d) contains the same number of elements as C{d)/l{d). The statement is 
now immediate. □ 

5. The Clifford Trace 

We now define the Clifford trace. The significance of this function for us is that 
every GP fiducial vector which has been constructed to date is an eigenvector of a 
Clifford unitary having Clifford trace = — 1. 

Let [F,x] G EC(d)/l{d) be the image of {F,x) under the homomorphism /e 
defined in Theorem |21 We refer to [F, x] as an extended Clifford operation (or 
Clifford operation if it S C(d) / 1(d)). The operators S [F, x] only differ by a phase. 
It is therefore convenient to adopt a terminology which blurs the distinction between 
the operation [F,x] and the operators U E [F,x]- In particular, we will adopt the 
convention that properties which hold for each U G [F, x] niay also be attributed 
to [-F, x]. Thus, we will say that [F,x] is unitary (respectively anti-unitary) if the 
operators U G [F,x\ are unitary (respectively anti- unitary). Similarly, we will say 
that 1-0) G is an eigenvector of [F, x] if it is an eigenvector of the operators 
U€[F,x]. 

It is easily verified that Tr(Fi) — Tr(F2) (mod d) whenever = [F'2,X2] 

(note that it is not necessarily true that Tr(Fi) = Tr(F2) (mod d) if d is even). 
We therefore obtain a well-defined function EC(d)/I(d) if we assign to each 

operation [F,x] the value Tr(F) (mod d). We obtain a function EC((i) Z^ by 
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assigning to each JJ G [F,x] the value Tr(F) (mod d). We use the term "Chfford 
trace" to refer to either of these functions. 

We now prove the main resuh of this section, which states that there is a con- 
nection between the order of a Chfford operation and its Chfford trace. 

Lemma 7. Let [F,x\ € C{d)/ l{d), where d is any dimension ^ 3. Then [F,x] 
of order 3 i/Tr(i^) = — 1 (modd). 

Let [F,x\ G where d is any prime dimension ^ 3. Then the stronger 

statement is true: [F,x] of order 3 if and only i/Tr(i^) — —1 {modd). 

Remark. The restriction to operations e C{d)/l{d) is essential (because if [F,x] is 
anti- unitary its order must be even). 

Proof. Let [F,x] G C{d)/l{d), and let k — Tr(F). Then, taking into account the 
fact that Det(F) — 1 (mod d), it is straightforward to show 

F^ ^ kF~1 (mod d) (115) 

implying 

F^ = {k^-1)F-k (modd) (116) 

1 + F + F^ = {k+1)F (modd) (117) 

Now suppose that k = — 1 (mod d). Then there are three possibilities: (a) d is odd; 
(b) d is even and k = — 1 (mod d); (c) d is even and k = —1 + d (mod d). In case 
(a) or (b) we have 

1 (modd) (118) 

1 + F + F^ = (modd) (119) 

while in case (c) we have — 1 = d^ — 2d ^ (mod d) , and consequently 

F'=(^V (modd) (120) 





1 + F + = (mod d) (121) 

Referring to the definition of Kj (see Theorem^ we deduce that, in every case, 

{F,xf = {FMl + F + F^)x)eL<f (122) 

implying that [F, x]'^ = [1,0]. It remains to show that neither [F,x] nor [F, x]^ = 
[1,0]. To see that [F,x] [IjO] observe that the contrary would imply — 1 = 
K = Tr(l) = 2 (mod d), which is not possible given that d ^ 3. Similarly, if 
[F, x]^ = [1,0] it would follow (taking the trace on both sides of Eq. (|115|l l that 
2 = K^ — 2 = — 1 (mod d), contrary to the assumption that d 7^ 3. We conclude 
that [F, x] is of order 3, as claimed. 

To prove the second part of the lemma suppose that d is a prime number 7^ 3 
and [F, x] is of order 3. Then (F^, {1 + F + F^)x) e Kf, implying F^ = 1 (mod d). 
In view of Eq. H116(l this means 

(k + 1) ((k- 1)F- 1) = (modd) (123) 

We now proceed by reductio ad absurdum. Suppose that k ^ —1 (mod d). Then 
Eq. H123|) and the fact that d is prime implies 

(k - 1)F = 1 (mod d) (124) 

Taking the trace on both sides gives (k + 1)(k — 2) = (mod d) implying k = 
2 (mod d). Substituting this value into Eq. (|124|l we deduce F = 1 (mod d), 
implying = 1 (mod d) and F^ ^ F (mod d). So 

{F,3x) = {F,xf eKf (125) 
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implying {F,x) S Kf. But that would mean [F,x] is of order 1, contrary to 
assumption. We conclude that k — —1 (mod d), as claimed. □ 

The result does not hold when d = 3 because then the identity has Clifford trace 
= — 1. It is, however, easily verified that in dimension 3 (as in every other prime 
dimension) every order 3 Clifford operation has Clifford trace = — 1. 

If d is not a prime number there may exist order 3 Clifford operations for which 
the Clifford trace 7^—1. Consider, for example. 



5 4 
2 -3 



-4 
5 



eC(6)/l(6) 



(126) 



Then [F, x] is of order 3 yet Tr(F) = 2 (mod 6). 

Because these results will play an important role in the following it is convenient 
to introduce some terminology. We will say that an operation [F,x] G C{d)/l{d) 
is a canonical order 3 unitary if 

(a) Tr(F) = -1 (mod d). 

(b) F is not the identity matrix. 

Note that the second stipulation is only needed because of the possibility that d = 3. 
li d ^ 3 an operation [F, x] G G{d)/ l{d) is a canonical order 3 unitary if and only 
if Tr(F) = -1 (mod d). 



6. The RBSC Vectors 

For 5 < d < 45 RBSC |17l have constructed CP fiducial vectors numerically. 
In this section we examine the behaviour of these vectors under the action of the 
extended Clifford group. In particular we show that each of them is an eigenvector 
of a canonical order 3 Clifford unitary. This suggests 

Conjecture A: CP fiducial vectors exist in every finite dimension. 
Furthermore, every such vector is an eigenvector of a canonical 
order 3 unitary. 

Conjecture A is related to a conjecture of Zauner's. Let 



[Z,0] 



(127) 



It will be observed that [Z,0] is defined, and £ C{d)/ l{d), for every dimension d, 
and that it is canonical order 3. Zauner has conjectured 

Conjecture B: In each dimension d there exists a GP fiducial 
vector which is an eigenvector of [Z, 0] . 

In Section [3 we will see that RBSC's numerical data also provides further support 
for Conjecture B. 

Let \^pd) be the RBSC vector in dimension d. In Table^we list, for each value of 
d, a unitary Clifford operation [Fd, Xd] having \tpd} as one of its eigenvectors. It will 
be seen that, in every case, Tr{Fd) — —1 (mod d), implying that [Fd, Xd] is canonical 
order 3. Clearly, ji/irf) is also an eigenvector of [Fd,Xd]'^- Moreover, also 
has Clifford trace = —1. There are, however, no other Clifford operations with 
these properties. 

In Table m we also list (n^i, ^^27 ^T-ds), the dimensions of the three eigenspaces of 
[Fd,Xd]i ^i^d Ud, the dimension of the particular eigenspace to which \ipd) belongs. 
It will be seen that, with one exception, \ipd) always belongs to an eigenspace of 
highest dimension (the exception being d = 17, where {ipd) belongs to the eigenspace 
of lowest dimension). 
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Xd {ndi,nd2,nd3) nd 



Fd 



Xd 



(ndi,nd2,nd3) "d 



12 



13 



14 



15 



16 



17 



19 



23 



24 



25 



3 



-2 
1 



-2 
1 



-10 
10 



(?) 



(9 



(-.') 



1,2,2 
1,2,3 
2, 2, 3 
2,3,3 
2,3,4 
3,3,4 
3,4,4 
3,4,5 
4,4,5; 
4,5,5 
4,5,6 
5,5,6 
5,6,6 
5,6,7 
6,6,7 
6,7,7 
6,7,8 
7,7,8 



7,8,9 



26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 



42 



43 



45 



-10 
-10 



-11 

-15 



7 

-10 



-31\ 
10 



-7 -; 

2 6 
12 



-29 
7 



-20 
21 



-3 
-12 



-10 
-6 



i-^ -/) (-3) 



-14 
6 



A) (ri) 



/-13 






I 16 




(■a 













(?) 








(-.'/ 


i) 










(-.' 




(-?o) 








(-0' 













(-.3 












(-A 


-;°) 
























y 





(8,9,9) 
(8,9, 10) 
(9,9, 10) 
(9, 10, 10) 
(9, 10, 11) 
(10,10,11) 
(10, 11, 11) 
(10, 11, 12) 
(11, 11, 12) 
(11,12,12) 
(11, 12, 13) 
(12, 12, 13) 
(12, 13, 13) 
(12, 13, 14) 
(13, 13, 14) 
(13, 14, 14) 
(13, 14, 15) 
(14, 14, 15) 
(14, 15, 15) 
(14, 15, 16) 



10 



10 



10 



12 



12 



12 



13 



13 



13 



15 



15 



15 



16 



Table 1. For each d the RBSC vector \tpd) is an eigenvector of 
the unitary operation Note that in every case TrF^ = 

— 1, implying that [Fd,Xd] canonical order 3. {ndi,nd2,nd3) 
are the dimensions of the three eigenspaces of [-Fd,Xd]) ^.nd na is 
the dimension of the eigenspace to which belongs. Note that 
nd = max(ndi, nd2, nas), with the single exception oi d= 17. 
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We used a computer algebra package (Mathematica) to construct the table. To 
illustrate the method employed we give a detailed description for the case d = 5. 
We begin with the observation that, if j-i/'s) is an eigenvector of [F, x], then 

(V'sI^pI^s) = e"^<>^^^P>(^5|^i=^p|^5) (128) 

for all p. So, using the value of \4>5) which is available on RBSC's website [53], we 
look for values of p, q such that 

|- (arg ((V-sI^pI^s)) - arg ((V'sl^qlV's))) (129) 

is an (approximate) integer. We find that if p = (1,0) this is only true when 
q = (1,0), (—1,1) or (0,-1) (mod 5), and that if p = (0,1) it is only true when 
q = (0, 1), (—1, 0) or (1, —1) (mod 5). Taking account of the requirement Det(F) = 
1 (mod 5) we deduce that the only candidates are (apart from the identity) 



1 -1\ 2 

1 r\2 



(130) 



and its square, [-FsjXs]^. To check that lips) actually is an eigenvector of [-F5,X5] 
we observe that is a prime matrix. So in view of Lemma0]we have the following 
explicit formula for the U G [F5, Xsl- 

U = ^e^'^(2,2) E e-'^^(^+^'-)|e.)(e.|^ (131) 

e'^ being an arbitrary phase. Suppose we choose 9 = j^. Then we find = 1 and 

||(J7-l)|7/^5>|r = (132) 

to machine precision. This confirms that j'f/'s) is indeed an eigenvector of [-FsjXsJ. 
To calculate the dimensions of the eigenspaces define, for r — 0, ±1 (and with the 
same choice oi 6), 

A- = ^ (l + e-^CZ + e^J/^) (133) 
Then Pr projects onto the eigenspace of (j with eigenvalue e^^. We find 

. - s fl r = l , , 

Tr(P,.)= „ ^ ^ (134) 

12 r = — 1 or 

implying that the dimensions of the eigenspaces are 1, 2, 2, and that \Tp^) is in one 
of the eigenspaces with dimension 2. 

In dimensions 6 to 45 the calculation goes through in essentially the same way. 
The calculation is, however, slightly more complicated when d is even, due to the 
fact that we must then require Det_Frf — 1 (mod 2d). Note, also, that when d = 
6, 21, 24, 28 or 36 the matrix Fd is non-prime, so we have to use the decomposition 
of Lemma |3| 

This method also enables us to establish the full stability group of \4>d)- i-e. the 
set of all operations (unitary or anti-unitary) e EC(c?)/l(d) of which \^d) is an 
eigenvector. It turns out that, with one exception, the stability group is the order 
3 cyclic subgroup generated by [Fd,Xd\- The exception is dimension 7, where 
the stability group is the order 6 cyclic subgroup generated by the anti-unitary 
operation 

iAr.(,]-[{\ -'),(;; 

Note that [^7,^7]2 = [F^^Xj]- 



(135) 



18 



d 


L 




Vd 


d 


L 


d 


Vd 


d 


L 


d 


Vd 


5 


(1 


OA 




19 


(" 


1 A 


(-•) 


33 




2 A 

-1.,) 












20 








34 








6 


(; 




(-.) 


(°. 




(-3) 




-n) 


(■3) 










21 








35 


/14 






7 




-3) 


& 


(i 




(:?) 


1,10 




(3) 






1 \ 


f-A 


22 


/I 






36 


fn 


1 \ 


(-A 


8 


(-1 


-3) 


(3) 




J 




(5 


-4) 


(-5) 


9 


(-3 




(") 


23 






(-^) 


37 


(-^. 




(=e) 


10 


( 3 


1 A 


fA 


24 


f 
(-1 


1 A 
-1) 


/'3\ 


38 


(-1 


1 A 

-J 


(isj 


11 


/I 


A 




25 


/I 


OA 


/'3\ 


39 




2A 


(n\ 




») 






J 


(0 


(2 




[u] 


12 


(-. 


i) 


& 


26 


(n 


-23^ 


(-) 


40 


(27 
1 14 


-35) 




13 








27 








41 








it 


») 






-0 




(-3 




(-■3) 


14 




1 \ 




28 


/ 12 


l\ 




42 


/ 2 


1 A 


^8^ 


(-. 


-1) 




1^-25 


26) 


(-3) 


(■. 


-36^ 


(') 


15 


(= 


-.) 


(?) 


29 


(" 




(-) 


43 


(-« 


-y 


&) 


16 




-■.) 




30 


(10 


3) 


(?) 


44 


1,-37 


20^ 


(i) 


17 








31 








45 


f 1 








3) 




{" 




(t) 


1,20 


;) 


(") 


18 




-«) 


(^) 


32 


( 27 

1-8 


■3) 


(-« 











Table 2. 



7. Zauner's Conjecture 

In the last section we saw that RBSC's numerical results support Conjecture 
A. Their results also support Conjecture B: i.e. Zauner's conjecture, that in each 
dimension d there exists a CP fiducial vector which is an eigenvector of [Z, 0] . 

In fact, for each 5 < d < 45 let [Ld,ri^] be the operation specified in Table [3 It 
is easily verified that 

[Ld,Vd][Fd,Xd][Ld,Vdr' = [Z,0] (136) 

This means that if [/ S [Ld,rij], and if \ipd) is the RBSC vector in dimension d, 
then Ultpd) is a CP fiducial vector which is an eigenvector of [Z, 0]. Conjecture B 
is thus confirmed numerically for every dimension < 45. 
This suggests 

Conjecture C: GP fiducial vectors exist in every finite dimension. 
Furthermore, every such vector is an eigenvector of a canonical 
order 3 unitary which is conjugate to [^,0]. 

Conjecture C is clearly stronger than Conjecture B. It also implies Conjecture A. 

An operation conjugate to [Z, 0] is automatically a canonical order 3 unitary. It 
would be interesting to know whether the converse is also true: i.e. whether every 
canonical order 3 unitary is conjugate to [Z, 0]. If that were not the case Conjecture 
C would be strictly stronger than Conjecture A. 
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Stability Group 




dimension 


type 


order 


number of orbits 


2 


non-Abelian 


6 


1 




non-Abelian 


6 


oo 


3 


non-Abelian 


12 


1 




non-Abelian 


48 


1 


4 


cyclic 


6 


1 


5 


cyclic 


3 


1 


6 


cyclic 


3 


1 


7 


cyclic 


3 


1 


cyclic 


6 


1 



Table 3. Stability groups in dimensions 2-7. In every case the 
stability group includes an order 3 cyclic subgroup generated by a 
unitary operation having Clifford trace = — 1. 



Dimensions 2 to 7: Vectors, Orbits and Stability Groups 



In dimensions 2-7 RBSC made a numerical search, in an attempt to find the 
total number of GP fiducial vectors. On the assumption that their search was 
exhaustive we use their data to calculate, for dimensions 2-7, the number of distinct 
orbits under the action of the extended Clifford group. We also calculate the order 
of the stability group corresponding to each orbit. Our results are tabulated in 
Table 13 They confirm that in dimensions 2-7 every GP fiducial vector is an 
eigenvector of a canonical order 3 Clifford unitary (in agreement with Conjecture 
A). We incidentally give exact expressions for two of the GP fiducial vectors in 
dimension 7 (one on each of the two distinct orbits). 

The calculations on which these statements are based are somewhat lengthy, and 
there is not the space to reproduce them here. We therefore confine ourselves to 
summarizing the end results, which it is straightforward (albeit tedious) to confirm 
with the help of (for example) Mathematica. 

Dimension 2. Exact solutions in dimension 2 have been obtained by Zauner [TB] 
and RBSC In dimension 2 the GP fiducial vectors all lie on a single orbit of 
the extended Clifford group. Consider the GP fiducial vector 



m = v^(3 + %/3)/6 |eo) + eT^(3-x/3)/6 |ei) 



(137) 



The stability group of 1-02) is the order 6, non-Abelian subgroup of EC(2)/I(2) 
generated by the unitary operation 



[F2,0] = 

and the three anti-unitary operations 

[^2,0] = 



1 

-1 -1 



1 

1 



(138) 



(139) 
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[B2,0] 



[C2,0] 







1 



1 

-1 -1 



(140) 



(141) 



Note that [^2,0] is canonical order 3. It follows from Lemmas El and El that 
|EC(2)/I(2)| 48. So the orbit consists of 48 ^ 6 = 8 fiducial vectors (identi- 
fying vectors wliich only differ by a phase), constituting 2 distinct SlC-POVM's (as 
described by RBSC). 

Dimension 3. Exact solutions in dimension 3 have been obtained by Zauner [TCj 
and RBSC (IJ- We saw in Section El that dimension 3 is unusual in that it is the 
only dimension for which the identity operator has Clifford trace = — 1. It seems 
to be unusual in another respect also: for it is the only case presently known where 
the CP fiducial vectors constitute infinitely many distinct orbits of the extended 
Clifford group. 

Consider the one parameter family of GP fiducial vectors 



1 



|e2)) 



(142) 



The complete set of GP fiducial vectors is obtained by acting on the vectors {ipsit)) 
with elements of EC(3). 

Let T and J be the operators defined by Eqs. © and (|1U6|I respectively. Then 



T|^3(t)) =-|^3(t + f)> 



and 



(143) 



So \ip3{t)) and \tp3{t')) are on the same orbit ii t' = ^ ±t for some integer n. At 
the cost of rather more computational effort one can show that this condition is not 
only sufficient but also necessary for \tp3{t)) and |V'3(*')) to the same orbit. 

So for each distinct orbit there is exactly one value of t G [0, ^] such that |V'3(0) is 
on the orbit. 

There are three kinds of orbit: a set of infinitely many generic orbits correspond- 
ing to values of t in the interior of the interval [0, and two exceptional orbits 
corresponding to the two end points t = and ^. 

The stability group of the exceptional vector |V'3(0)) consists of all 48 operations 
of the form [F, 0], where F is any element of ESL(2, Z3). The orbit thus consists 
of 432 48 = 9 fiducial vectors, constituting a single SIC-POVM. 

The stability group of the exceptional vector |'!/'3(f)) is the order 12 non-Abelian 
subgroup of EC(3)/I(3) generated by the unitary operation 



[^^3,X3] = 

and the anti-unitary operation 

[^3,X3] = 

Note that 

[^3,^3]' = 



-1 
-1 







1 

-1 



-1 
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(145) 



12 



(146) 

36 fiducial vectors, 



is canonical order 3. The orbit thus consists of 432 
constituting 4 distinct SIC-POVMs. 

The stability group of a generic vector \ip3{t)) with < < < ^ is the order 6 
non-Abelian subgroup generated by the unitary operation 



[F3,X3? 



1 

-1 



(147) 
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and the anti-unitary operation 



-1 

-1 1 



(148) 



The orbit thus consists of 432 6 = 72 fiducial vectors, constituting 8 distinct 
SIC-POVMs. 



Dimension 4. The vector 



75 



40 



2 cos 



J|eo) +^(e-'^ + (2 + Vs) ^e^) |ei) 
+2isin-|e2) + i{e^'^ - (2 + \/5)' 



(149) 



is a GP fiducial vector in dimension 4, as discovered by Zauner flF and RBSC |17p. 
The stabihty group of IV'4) is the order 6 cyclic subgroup of EC(4)/I(4) generated 
by the anti-unitary operation 



[A4, Xi] = 



Note that 



1 

-1 3 



-1 1 
-1 2 



(152) 



(153) 



is canonical order 3 (where we used Eq. (|81|l to obtain the last expression on the 
right hand side). 

It follows from Lemmas El and El that the group EC(4)/I(4) is of order 1536. So 
the orbit generated by \tlj4) contains 1536 6 = 256 fiducial vectors, constituting 
256 16 = 16 SIC-POVMs. It was shown by RBSC that there are only 16 SIC- 
POVMs in dimension 4. We conclude that the fiducial vectors all lie on a single 
orbit of the extended Clifford group. 



Dimension 5. Let jf/'s) be RBSC's numerical vector in dimension 5. We noted 
in the last section that the stability group of jV's) is of order 3. It follows from 
Lemmas Eland El that the group EC(5)/ 1(5) is of order 6000. So the orbit generated 
by li/^s) contains 6000 3 = 2000 fiducial vectors, constituting 2000 25 = 80 SIC- 
POVMs. It was shown by RBSC that there are only 80 SIC-POVMs in dimension 5. 
We conclude that the fiducial vectors all lie on a single orbit of the extended Clifford 
group. 

Note that Zauner's analytic solution in dimension 5 (on p. 63 of his thesis jl6| 'l 
can be used to give exact expressions for each of the vectors on the orbit. 



In Zauner's notation I1/14) is the vector 

e-t {XtPia + P^YtPu) (150) 

In RBSC's notation it is the vector 

ro|eo> + r+e'*+ |ei> + ne'"^ |e2> + r.e**- |e3> (151) 
for the case n = j = m = 1 and fc = (note, however, that there is a typographical error in 
RBSC E): their expression for ro should read ro = W (l - l/Vf>) / 2 - V2) ). 
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Dimension 6. Let jV'e) be RBSC's numerical vector in dimension 6. We noted 
in the last section that the stability group of jV'e) is of order 3. It follows from 
LemmasElandElthat the group EC(6)/ 1(6) is of order 10368. So the orbit generated 
by l^/ig) contains 10368-;- 3 — 3456 fiducial vectors, constituting 3456-^-36 = 96 SIC- 
POVMs. It was shown by RBSC that there are only 96 SIC-POVMs in dimension 5. 
We conclude that the fiducial vectors all lie on a single orbit of the extended Clifford 
group (in agreement with Grassl's j21| analysis, based on his exact solution in 
dimension 6). 

Note that Grassl's Jl, analytic solution can be used to give exact expressions 
for each of the vectors on the orbit. 

Dimension 7. Let be RBSC's numerical vector in dimension 7. We noted 
in the last section that the stability group of jV'y) is of order 6. It follows from 
LemmasElandElthat the group EC(6) / 1(6) is of order 32928. So the orbit generated 
by IV't) contains 32928 6 = 5488 fiducial vectors, constituting 5488-^49 = 112 
SIC-POVMs. However, it was shown by RBSC that there are 336 SIC-POVMs in 
dimension 7. We conclude that there must be at least one other orbit. 

The search for the additional orbit or orbits is facilitated by the fact that in 
dimension 7 there exists a canonical order 3 Clifford unitary for which the F matrix 
is diagonal: namely 



The fact that Flj is diagonal means that the IJ G [Fy, 0] are permutation matrices. 
Specifically 



for every [/ e \F^^ 0] (where e** is an arbitrary phase, and where we have used the 
decomposition described in LemmaOJ. This considerably simplifies the calculations. 
We will also have occasion to consider the anti-unitary operation 




(154) 



6 




(155) 




(156) 



which is a square root of \Flj , 0] . 

We look for eigenvectors of 0]. Let 




if r 



if r 



1,2 or 4 
3, 5 or 6 



(157) 



Also let 




(158) 



and 




(159) 
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Then define 



1^7) = ao|eo) 



(ai + lra2) \er) 

r=l 



m=bo\eo)+^b 



le 



(160) 



(161) 



It is readily confirmed that Itp'^) and {ipj) are both GP fiducial vectors. The stability 
group of 1 7/17) is the order 3 subgroup generated by [F^, 0], while the stability group 
of li/)") is the order 6 subgroup generated by [Ay,0]. Since the stability groups 
are non-isomorphic the orbits generated by {ip'-^) and Itp") are disjoint. The orbit 
generated by \il>'^) contains 32928^3 = 10976 fiducial vectors, constituting 10976-^ 
49 = 224 SIC-POVMs. The orbit generated by \ip'-l) contains 5488 fiducial vectors, 
constituting a further 112 SIC-POVMs. This accounts for aU 336 of the SIC-POVMs 
identified by RBSC. We conclude that there are no other orbits, apart from these 
two. 

For the sake of completeness let us note that 



1^7) = W) 



where {ipj) is RBSC's numerical vector and ?7 is a unitary operator 

U e 



1 1 

-3 -2 



(162) 



(163) 



Finally, let us remark that Ir is the Legendre symbol (see, e.g., Nathanson 
or Rose j34j) 

It has the important property that Irs — Irh for all r, s G Z. 



(164) 



9. A Fiducial Vector in Dimension 19 

In Section we saw that, except in dimension 7, each of RBSC's numerical 
solutions has stability group of order 3. This might encourage one to speculate 
that when d > 7 the stability group is always of order 3. In this section we show 
that there is at least one exception to that putative rule, by constructing a GP 
fiducial vector in dimension 19 for which the stability group has order > 18. 

The vector we construct is an eigenvector of the order 18 anti-unitary operation 



[^'19, 0] = 



-9 \ /G 
-2j'U 



Note that 



e EC(19)/I(19) 



-8 0\ fO 
7/ ' lo 



(165) 
(166) 



is canonical order 3. 

The construction is similar to our construction of the vector j?/'") in the last 
section. Let l'^ be the Legendre symbol 



1 



if r = 1,4,5,6,7,9, 11, 16 or 17 
if r = 2, 3, 8, 10, 12, 13, 14, 15 or 18 



and let 



b'o = 



'5 + 9V5 /lO-x/5 ^, _i 

61 = 1/ — —. — = cos 



'V5-1 



95 



190 



(167) 



(168) 
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Then define 

18 

IV'i9)=6o|eo)+E^i^''^''l^'-) (169) 

r=l 

It is readily confirmed that IV'ig) is a GP fiducial vector, and an eigenvector of 

[^'l9,0]. 

Observe that the orbit generated by iV'lg) is disjoint from the orbit generated by 
RBSC's numerical vector iV'ig) (because the stability groups are non- isomorphic) . 
It follows that there are at least two distinct orbits in dimension 19. 

10. DiAGONALIZING THE F MATRIX 

Our construction of the exact solutions \'ip'-j), IV'?) and ji/'ig) in Eqs. H16U|) . (|161|l 
and (|169|l was facilitated by the fact that in dimensions 7 and 19 there exist canoni- 
cal order 3 unitaries for which the corresponding F matrix is diagonal. It is natural 
to ask in what other dimensions that is true. The theorem proved below answers 
that question. 

We will need the following lemma: 

Lemma 8. Let p be a prime number = 1 {mod 3), and let n be any integer > 1. 
Then there exists an integer a such that 

a^ + a + l = {modp'^) (170) 

Proof. The proof relies heavily on the theory of primitive roots, as described in (for 
example) Chapter 3 of Nathanson or Chapter 5 of Rose Let (/) be Euler's 
phi, or totient function (so for every integer a: > 1, <j){x) is the number of integers 
y in the range 1 < y < x which are relatively prime to a;). Then there exists a 
single positive integer g such that for every integer m > 1 the multiplicative order 
of g, considered as an element of Zp™, is 0(p™) = {p — 1)]3'"~^ (see, for example, 
Nathanson P- 93, or Rose "37, p. 91). The fact that p = 1 (mod 3) means 
p = 3A; + 1 for some integer A; > 1. Define 

a = /f"" (171) 

It is then immediate that 

It is also true that a — 1 is relatively prime to p. For suppose that were not the 
case. It would then follow from the definition of a, and the fact that g is a primitve 
root modulo p, that 

kp''-^ ^ l{p - \) ^ m (173) 

for some integer Z > 1. That, however, is impossible since p is not a multiple of 3. 

The fact that a — 1 is relatively prime to p means that there exists an integer (3 
such that 

(i{a - 1) = 1 (mod p") (174) 
It now follows from Eqs. ifTT^ and (|T7^ that 

0,2 +a+ 1 = /3(a3 - 1) = (mod p") (175) 

□ 

We are now in a position to prove our main result: 

Theorem 3. There exists a canonical order Z unitary [F,'x\ G C{d)/l{d) for which 
the matrix F is diagonal if and only if each of the following is true 

(1) d has at least one prime divisor — 1 {mod 3). 

(2) d has no prime divisors — 2 {mod 3). 

(3) d is not divisible by 9. 
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Remark. So there exist canonical order 3 unitaries [F, x] for which F is diagonal in 
dimension 7, 13, 19, 21, 31, 37, 39, 43, 49, . . . 

Proof. We begin by proving sufficiency. Suppose that conditions (1), (2) and (3) 
are all true. Then we have, for some t > 1, 

d = r"pi'...pT (176) 

where the Pi are distinct prime numbers = 1 (mod3), where the integer no = or 
1, and where the integers ni, . . . , nt are all > 1. It follows from Lemma|Hlthat there 
exist integers ai,. . . ,at such that 

a- + a, + 1 = (modp;"') (177) 

for i — l,...,t. We then use the Chinese remainder theorem (see, for example, 
Nathanson or Rose Pl]) to deduce that there exists a single integer a such 
that 



and 



for z = 1, . . . , t. We have 



and 



a = 1 (mod 3) (178) 

a = ai (modp"') (179) 
+ a + 1 = (mod 3) (180) 



a2+a + l = (modp"') (181) 
for i = 1, . . . , t. Consequently 

+ a + 1 = (mod d) (182) 

It follows that the matrix 

e SL(2, Z^) (bearing in mind that d is odd). Moreover, Tr(i^) = —1 (mod d). Since 
d 7^ 3 we conclude that [F, x] is a canonical order 3 unitary for all x G {'^df' ■ This 
proves sufficiency. 

To prove necessity suppose 

F= (^J °) eSL(2,Z3) (184) 

is such that [F, x\ is canonical order 3 for some x ^ (^d)^- Then a+5 = — 1 (mod d), 
implying 

+ a + 1 = (mod d) (185) 
= 1 (mod d) (186) 

(in view of the fact that aS — \ (mod d)). 

To show that d has no prime divisors — 2 (mod 3) assume the contrary. It would 
then follow from Eqs. ifT^ and itTF^ that 

+ a + 1 = (mod p) (187) 

a^ = l (modp) (188) 

for some prime number p — 2 (mod 3) . Let r be a primitive root of p and let fc € Z 
be such that < < p — 1 and a = r'' (mod p) (see, for example, Nathanson 
or Rose |311)- Then Eq. (|188|l implies r^'^ = 1 (mod p) which, in view of the fact 
that r is a primitive root, means 3fc ~ l{p — 1) for some / G Z. The fact that 
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< k < p — 1 implies < I < 3. Taking into account the fact that p — 1 is not 
divisible by 3 (because p = 2 (mod 3)) we deduce that / — 0. But then fc = 0, 
implying a = 1 (mod p). In view of Eq. H187(l this means 3 = (mod p): which is 
a contradiction. 

To prove that d is not divisible by 9 we again proceed by reductio ad absurdum. 
Suppose that d were divisible by 9. It would then follow from Eq. H185|) that 

+ a + l = {) (mod 9) (189) 

However, it is easily verified (by explicit enumeration) that this equation has no 
solutions. 

Finally, suppose that d had no prime divisors = 1 (mod 3). In view of the 
results just proved it would follow that d — i. But if = 3, Eq. (|185|l implies 
a = 1 (mod 3). Taking into account the requirement a5 = deti^ = 1 (mod 3) this 
means i5 = 1 (mod 3). But then F is the identity matrix, which contradicts the 
assumption that ~x\ is a canonical order 3 unitary. We conclude that d must 
have at least one prime divisor = 1 (mod 3). □ 

11. Conclusion 

RBSC conclude their paper by saying "a rigorous proof of existence of SIC- 
POVMs in all finite dimensions seems tantalizingly close, yet remains somehow 
distant" . That well expresses our own perception of the matter. While working on 
this problem we have several times had the sense that the crucial discovery lay just 
round the corner, only to find that our hopes were illusory. We make our results 
public in the hope that they may, nevertheless, contain a few clues, which will help 
to take us further forward. 

In particular it seems to us that significant progress would be made if it could be 
established whether it is in fact true that every GP fiducial vector is an eigenvector 
of a canonical order 3 unitary. Also, if that is the case, one would like to know 
exactly why it is the case. 
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